Gelling by Heating 



m 
o 



O 



o 
o 



> 

(N 
00 

in 
(N 

cn 
o 
m 



x 



Sandalo Roldan- Vargas 1 , Frank Smallenburg 1 , Walter Kob 2 , Francesco Sciortino 1 
1 Department of Physics, Sapienza, Universita di Roma, Piazzale Aldo Moro 2, 1-00185, Roma, 
2 Laboratoire Charles Coulomb, UMR 5221, CNRS and Universite Montpellier 2, Montpellier, France 

We introduce a simple model, a binary mixture of patchy particles, which has been designed to 
form a gel upon heating. Due to the specific nature of the particle interactions, notably the number 
and geometry of the patches as well as their interaction energies, the system is a fluid both at high 
and at low temperatures T, whereas at intermediate Ts the system forms a solid-like disordered 
open network structure, i.e. a gel. Using molecular dynamics we investigate the static and dynamic 
properties of this system. 

PACS numbers: 83.80.Kn, 47.57.-s, 64.70.qd, 61.20.Ja 



Some of the most versatile and efficient strategies for 
designing new materials with unconventional behavior 
are based on the idea of competitive interaction. While 
in biology the term competitive interaction indicates the 
result of rivalry between two or more species compet- 
ing for resources, in physics it stands for the presence of 
several interaction mechanisms that can stabilize com- 
peting local structures, leading to novel and highly in- 
teresting features of the system. Recent examples for 
this mechanism in soft matter systems include, among 
others, the competition between short-range attraction 
and long-range repulsion in charged colloids giving rise 
to cluster phases PHI , the competition between chaining 
and branching in patchy colloids 0, [E[ where a specific 
design of the inter-patch interactions results in a phase 
diagram in which the density of the coexisting liquid ap- 
proaches the density of the gas and the design of 
DNA-coated colloids with two different DNA sequences 
for the purpose of establishing a competition between 
intra and inter-particle interactions, favoring crystal for- 
mation Q. Often the very nature of these competing 
mechanisms promotes the emergence of a structure con- 
trolled by energy (stable at low T) which competes with 
a structure stabilized by entropy at intermediate T. 

Recent progress in the synthesis of colloids have led to 
a new generation of particles with highly directional and 
selective interactions, providing valence to colloids p"}- 
The ability to tune the interactions between nano- and 
meso-sized particles almost at will [8|, |a|, a nd to design 
the geometric properties of the patches [I0| and/or their 
functionalization, offers today the possibility to exploit 
the idea of competitive interactions to modulate mate- 
rial properties with external control parameters. In this 
Letter we develop one such possibility: The design of a 
material whose viscosity increases upon heating. In par- 
ticular, we show how a simple design of a binary mixture 
of limited valence particles can indeed provide a model 
where the competition between entropy and potential en- 
ergy causes the system to show a re-entrant behavior, 
passing reversibly from a fluid to a gel and again to a 
fluid when T is varied. 

It has been shown that patchy colloidal particles with 
a limited number / of attractive patches progressively 



cluster when cooled, leading such system to form a perco- 
lating network that at sufficiently low T will incorporate 
all particles in the system [llT - [f6j . During this process, 
the diffusion coefficient progressively decreases and the 
viscosity simultaneously increases. In the network state 
the lifetime of the bonds (a T-controlled quantity) fixes 
the timescale over which the system behaves as a solid. 
It is known that when / = 4, the particles (here called A 
species) form a random tetrahedral network which closely 
resembles the structure of network-forming atomic sys- 
tems like silica and silicon [I?} • To melt the network at 
low T, we consider the addition of a second species (B) 
with a single bonding patch that competes for bonding 
with the patches on the network forming A-species. The 
idea is to design a competitive mechanism such that the 
bonding between A- and B-particles becomes dominant, 
but only at a temperature much lower than the one at 
which the ^4^4-network is formed. As a result, the stable 
low-T phase consists of A-particles decorated with / B- 
particles which are free to diffuse in the sample volume, 
whereas at intermediate T the system forms a highly vis- 
cous AA-network that is progressively fragmented and 
transformed into a fluid upon heating. 

Model: We consider a binary mixture of patchy col- 
loids where each A-particle has / = 4 patches on its sur- 
face that are arranged in a tetrahedral geometry and the 
B-particles have only one patch (see Fig. [TJ. The patch- 
patch interaction is modeled via a Kern-Frenkel poten- 
tial [Ilj], a model that has been extensively used over 
the last decade to compare simulations and experiments 
on the self-assembly of patchy colloids fl9i - f2Tl j . Each A- 
patch can interact either with a B-patch with unit energy 
cab and bonding volume Vab or with another A-patch 
with energy €aa — 0.956^ and Vaa (3> Vab)- No BB- 
bonding is allowed. The bonding volumes are each deter- 
mined by an interaction range 5 a and an angular patch 
width 6 a (a G {AA,AB}) (see Fig. [T] and Ref. 
The attractive patch-patch interaction is complemented 
by an isotropic hard-core repulsion, where the spherical 
cores have diameters a a and o~b — 0.35avi. The size ra- 
tio was chosen such that the B-particles can block the 
A-patches from bonding to other A-patches without sig- 
nificantly increasing the packing fraction of the pure A 



2 



a) b) 




c) d) 




FIG. 1. a) Schematic of the interaction parameters in the 
Kern-Frenkel model. An A-patch can bond with either an- 
other A-particle (b), or with a £>-particle (c). Panel (d) 
shows a "flower", i.e. a fully bonded cluster consisting of 
one A-particle and four B-particles, representing the lowest- 
energy state of the system. Here the interaction ranges and 
the angular patch widths are cosOaa = 0.92, Saa = 0.15ita, 
and cosOab = 0.99, Sab = 0.2<ta- With these choices 
Vaa = 3.49 • 10~ 3 cri and Vab = 3.79 • 10" 5 ai 0. Note 
that in the binary Kern-Frenkel model, 9 and S of each bond 
are defined by the species of both bonding partners and are 
not properties of individual particles. 



system. Due to geometric constraints, each patch can be 
involved in only a single bond. 

Using event-driven molecular dynamics simulations 
[22T - I25T ] . we have studied a system of Na — 600 and 
Nb = 2400 particles, corresponding to a total num- 
ber density po\ — 3.0, with partial number densities 
Pa&a = 0.6 and Pbo~ a ~ 2.4, for a wide range of T, 
whose unit is given by cab/^b, where ks is Boltzmann's 
constant. The composition of the system is thus fixed 
at xa — 0.2. The density pa&\ = 0.6 °f A-particles 
corresponds to the optimal density at which tetrahedral 
particles form an unstrained fully bonded network (2fjj. 
With this composition, the fully bonded network has an 
energy of 2Na^aa, whereas a configuration in which all 
the B-particles are bonded to the A-particles has a sig- 
nificantly lower energy of ANacab- We also simulate for 
a low T (= 0.0AeAB/kB) a reference system composed of 
600 flowers (see Fig. [T}1), i.e. A-particles bonded to four 
-B-particles. At this low T no bond breaking events take 
place within the simulation time. 

Results: Figure [2] demonstrates the basic mechanism 
of the competitive interactions present in our system. It 
shows the probability that a patch on an A-particle is 
bonded to another A-patch (paa) or to a i?-patch (pab) 
as a function of the inverse T. On cooling, paa starts 
to grow, signaling the onset of the network formation, 
reaching a maximum around T — O-lle^s/fce- We recall 
that within a mean-field description, percolation of par- 
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FIG. 2. Probability that a patch on an A-particle is bonded 
to another patch on an A-particle, paa (black circles), or to 
a patch on a B-particle, pab (red squares). The dashed lines 
are the prediction for these probabilities as obtained from the 
Wertheim theory. 



tides with valence / = 4 takes place at paa = 1/3 [13] ■ 
Since at the maximum we find paa ~ 0.9 we can con- 
clude that at this T the A-particles have formed a highly 
bonded percolating network. Upon further cooling, pab 
significantly increases, showing that the AB-bonds are 
starting to replace the AA ones, i.e. while pab ap- 
proaches 1.0, paa — > 0. The entropy associated with 
the larger bonding volume Vaa for the AA interaction 
is crucial for promoting the formation of a large num- 
ber of AA-bonds at intermediate T, before the energeti- 
cally preferred but entropically disfavored AB-bonds set 
in. Fig. [2] also shows the parameter-free theoretical pre- 
dictions for the bonding probabilities as obtained from 
the first-order thermodynamic perturbation theory de- 
veloped by Wertheim [28 M30J ( details on the Wertheim 
calculations are reported in [22|). The Wertheim theory 
nicely captures the mechanism of competing interactions, 
reproducing the position and height of the maximum of 
Paa as well as the low T trends of paa and pab • 

Figure [3] shows the unusual T dependence of the struc- 
ture of the system, which also reflects the non-monotonic 
behavior of paa- At high T, the partial structure fac- 
tor Saa{q) shows the conventional q— dependence found 
in simple liquids with a main peak around q<JA = 7.2. 
Upon decreasing T the main peak splits into two, one 
located around qo~A = 5.2 and a higher one around 8.4. 
This double peak feature is typical of liquids that have a 
local tetrahedral network structure, such as silicon or sil- 
ica [l7|. The peak at qoA ~ 8.4 corresponds to the near- 
est neighbor distance between two bonded A-particles, 
whereas the one placed around 5.2 is associated with 
the second-nearest neighbors in the tetrahedral network. 
Note that this double peak structure is most pronounced 
at T w O.lleAs/fcs, i.e. at the T at which paa has a 
maximum (see Fig. [5]) and hence the gel is maximally 
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FIG. 3. Partial structure factor Saa(?) for different values 
of T (solid lines with different symbols). The Saa{(}) for a 
fluid of flowers at T — O.OieAB /fes (black solid line) is also 
represented. 



connected. When T is lowered even further the dou- 
ble peak structure disappears and Saa (<?) becomes again 
similar to the structure factor of a fluid composed of flow- 
ers (which is also included in Fig. [3] as a reference). Since 
the size of a flower is larger than that of an A-particle, the 
peak position at low T is to the left of the one observed 
at high T. 

We now quantify the effect of the competing interac- 
tions on the dynamics of the system and provide evidence 
that the change of the structure rich in AB-bonds at low 
T to the highly bonded AA-network generates a slowing 
down of the dynamics on heating. To do this, we calcu- 
late the mean squared displacement (MSD) for the parti- 
cles of both species and then their corresponding diffusion 
coefficients D a (a G {A,B}) from the long-time behavior 
of the MSD via the Einstein relation [31(. To subtract 
the trivial trend originated from the T-dependence of the 
thermal velocity we divide D a by a reference diffusion co- 
efficient Dq = ct a /to, where tq = y/ thaCa/ ^bT and rriA 
is the mass of an A-particle. 

Figure |4] shows the T-dependence of D a in an Arrhe- 
nius plot. At high T, D a is approximately constant for 
both type of particles, indicating that bonds do not play 
a significant role. On cooling, Da starts to decrease very 
rapidly, with a super- Arrhenius T-dependence reminis- 
cent of that observed in molecular networks (l7j , turning 
into an Arrhenius law with an activation energy approx- 
imately equal to 2caa (see dashed line in Fig. HJ. Sim- 
ilar values of the activation energy are typically found 
in tetrahedral network-forming systems where most of 
the particles belong to the percolating cluster, and bond 
breaking is the bottleneck for relaxation [l5|, [H, |33| . Be- 
fore the gel starts to decompose at a temperature be- 
low T w O.lleAs/fcs, Da has already decreased by four 
orders of magnitude compared to its value at high T, 
indicating the formation of a persistent network. For 
T < O.lle^s/fcsj Da starts to increase. This rising per- 
sists down to the lowest T at which we were able to equili- 
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FIG. 4. Arrhenius plot of the normalized diffusion coefficient 
D a /Do (a — A, circles, a = B, squares). Also included is an 
Arrhenius law with activation energy 2eaa (dashed line) . The 
diffusion coefficient of a fluid of flowers at low T is represented 
by a blue diamond. 



brate the system. We emphasize that this non-monotonic 
T-dependence is only observed for the A-particles, i.e. 
the particles which are involved in the formation of the 
network. In contrast, Dg shows only a rather mild T- 
dependence. Figure H] also shows the diffusion coefficient 
of the fluid of flowers at T = 0.04eAs/fcs for which 
Da = Db- This common value of the diffusion coeffi- 
cient provides a low-T limit to which both Da and Db 
converge, consistent with the trends shown by Da and 
D B at low T. 

A non-monotonic behavior of the characteristic time 
is also found in the time evolution of the collective- and 
self-intermediate scattering functions. Their study pro- 
vides insight into how the relaxation dynamics depends 
on the considered length scale. Figure [5^ shows an Ar- 
rhenius plot of the relaxation time ta {q) determined from 
the time integral of the intermediate scattering function 
of the ^4-particles for two different q— vectors: q<JA = 5.2 
and qoA — 7.2, which correspond, respectively, to the lo- 
cation of the first peak in the network and in the high-T 
fluid (see Fig. [3]). We find that the self and collective 
relaxation times TA{q), normalized by To, show qualita- 
tively the same T-dependence: a plateau at high T, a 
fast increase within the T-range in which the network 
is formed, a quick decrease once the network starts to 
break up again, and a final plateau at low T. This T- 
dependence is observed for both values of q, indicating 
that the relaxation mechanism does not depend on the 
length scale considered. Analogous to the diffusion coef- 
ficient, the self and collective dynamics of the B-particles 
are found to be faster than those of the ^4-particles and 
their q and T dependence will be reported elsewhere. 

To provide further evidence that the system is ergodic 
on long time scales, i.e. that the structure of the system 
has completely lost its memory of the initial state, we 
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FIG. 5. a) Arrhenius plot of the normalized relax- 

ation time TA(q)/ro as obtained from the self (dashed lines) 
and collective (solid lines) scattering functions, where tq = 
^/m A cr 2 A /kBT. The different curves correspond to the wave- 
vectors given by the first two peaks in Saa{<})- b) Arrhenius 
plot of the normalized bond-persistence time it, /to for the 
A A (squares) and AB bonds (circles). Also included is an 
Arrhenius law with activation energy cab (blue solid line). 



investigate the bond persistence function p&(t), i.e. the 
probability that a bond which is present at time zero is 
also present at time t. When Pb(t) approaches zero, all 
bonds which were present at time zero have been bro- 
ken. Hence, the relaxation time of pb{t) provides infor- 
mation on the restructuring time of the network con- 
nectivity. Figure [5]d shows an Arrhenius plot with the 
T-dependence of the decay time Tf,, where Pb(rb) = e _1 . 
At intermediate and low T, Tb is larger than the relax- 
ation times shown in Fig. [5^- We find thus that Pb(t) 
decays to zero only on a time scale that is significantly 
longer than the relaxation times associated with the scat- 
tering functions, confirming that some fraction of spacial 
decorrelation of the network, as quantified by the col- 
lective scattering function, takes place at partially fixed 
bonding pattern. In other words, while the decay to zero 
of the scattering functions for a given q— value implies 
that the particles have moved over a distance on the or- 



der of 2n/q, such motion does, however, not necessarily 
require that all bonds are broken, since, e.g., a cluster of 
particles can move in a collective manner. A bond can 
persist up to very long times even if the structure of the 
system changes significantly. The extreme case occurs at 
very low T where the system is a fluid of flowers that 
relaxes relatively quickly but in which AB-bonAs survive 
for a very long time. Figure [SJj also shows the differ- 
ent T-dependence of the lifetime of the AA- and of the 
AB-bonds. An Arrhenius dependence with an activa- 
tion energy very close to cab (see dashed line in Fig. [SJd) 
characterizes the Ai?-bonds in the entire T-range, sug- 
gesting that the mechanism for the breaking of an AB- 
bond is not collective in nature but a simple activated 
process. In contrast to this behavior, the breaking time 
for an AA-bond follows an Arrhenius law at high T but 
becomes super- Arrhenius within the T-region in which 
the gel forms, showing that the bond-breaking process 
becomes coupled to the degree of bonding. At low T, 
when most of the A-sites are bonded to i?-particles and 
the network is disrupted, the effective activation energy 
decreases again to recover at very low T an Arrhenius 
behavior with an activation energy given by £aa- 

Conclusions: In summary, we have shown that a ju- 
dicious choice of the interaction parameters of a binary 
mixture of A and B patchy particles allows us to gener- 
ate a non-monotonic T-dependence of its dynamic prop- 
erties. Essentially, we set up a competition between 
network-forming AA-bonds and network-breaking AB- 
bonds, and tune the balance between both bond types 
by choosing their bonding volumes (and therefore the 
entropy) and energies. At high T, very few bonds are 
formed, and the system behaves similar to a binary fluid 
of hard spheres. At slightly lower T, the stability of a 
stiff, percolating network is ensured by the larger en- 
tropy associated with the much larger bonding volume 
Vaa- However, when T is decreased even further, the 
system instead forms a fluid of small clusters, stabilized 
by a much larger number of AB-bonds, corresponding 
to a lower potential energy. We have shown that, com- 
pared to both the high-T and low-T fluids, the network 
state relaxes and diffuses more slowly by several orders of 
magnitude. Thus, the system forms a reversible gel that 
transforms into a fluid upon both heating and cooling. 
We conclude by listing two experimental systems which 
we believe are very promising candidates to experimen- 
tally test the ideas presented in this Letter: a solution 
of DNA constructs of valence four [§3, [35| in the pres- 
ence of competing DNA single strands and the binary 
mixture of patchy particles recently synthesized Q. Both 
these model systems have the potential to provide soft 
materials that gel on heating. 
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